Intense terahertz laser fields on a two-dimensional hole gas with Rashba spin-orbit 

coupling 
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We investigate the influence on the density of states and the density of spin polarization for a 
two-dimensional hole gas with Rashba spin-orbit coupling under intense terahertz laser fields. Via 
Floquet theorem, we solve the time-dependent Schrodinger equation and calculate these densities. 
It is shown that a terahertz magnetic moment can be induced for low hole concentration. Different 
from the electron case, the induced magnetic moment is quite anisotropic due to the anisotropic 
spin-orbit coupling. Both the amplitude and the direction of the magnetic moment depend on the 
direction of the terahertz field. We further point out that for high hole concentration, the magnetic 
moment becomes very small due to the interference caused by the momentum dependence of the 
spin-orbit coupling. This effect also appears in two-dimensional electron systems. 
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Optical properties of semiconductors are sensitive to 
external conditions. Almost fifty years ago Franai and 
Keldysh 2 pointed out that under static electric fields the 
absorption coefficient becomes finite below the band gap, 
and the above-gap absorption spectrum shows oscilla- 
tions. In the late 1990s, Jauho and JohnsenM studied the 
optical properties of semiconductors under strong ac field 
and developed dynamic Franz-Keldysh effect (DFKE), 
which presents the blueshift of the main absorption edge 
and the fine structure near the band gap. This ef- 
fect is particularly obvious for semiconductors under in- 
tense terahertz (THz) field, thereby leading to extensive 
theoretical and experimental interests on THz electro- 



is small enough so that only the lowest subband is rel- 
evant. For this structure, the lowest subband is heavy 
hole (HH) like. By applying a suitable strain, it can be 
light hole (LH) like. The confinement is assumed to be 
large enough so that the lowest HH and LH subbands 
are well separated and we can consider the HHs and LHs 
separately. The Hamiltonian can be written as (h = 1) 



H x (K,t) = 



K 2 

2m* 



(1) 



Here A = LH, HH and K = k — eA(t), with k standing 



and can be written as if* (K) = ^(a x rt x +a y tty). Here, 
a is the Pauli matrix. For HHs, we hav o 13 i 14 
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K) = 2K X (^K 



(2) 
(3) 



optics£&I£ Very recently Cheng and Wu brought the for the electron momentum. m\ is the effective mass, 
spin degrees of freedom into the study of the THz field For GaAs Q w > H *o is mainly due to the Rashba term^ 
induced effects. They studied a two-dimensional electron 
gas (2DEG) with the Rahhba spin-orbit coupling (SOC) 
under intense terahertz fields It is shown that the THz 
field can efficiently modify the density of states (DOS) of 
the 2DEG and excite a magnetic moment oscillating at 
THz frequency. Later, Jiang et al. studied similar effects 
of quantum dots 1 ^, and further discussed the spin dis- 
sipation under THz driving fields^ However, all these 
works concentrate on electron systems. Up to now there 
is no study on the spin properties of hole systems under 
intense THz field. In this letter we study the effect of in- 
tense terahertz laser fields on a two-dimensional hole gas 
(2DHG) with Rashba SOC and show that this system 
has some new properties different from the previously 
studied electron system. 



with j« H = E z {^ h + 75 7 4 "'>), j?" 
3~fl% 7h ). For LHs, 



ni H (K) = 2K y {^ H Kl+^ H Kl+^ H {kl)) ,(4) 



n^(K)=-2K x ( 7 ^K' x 



7r<*s»,(5) 



We consider ap-type GaAs (100) quantum well (QW). 
The growth direction is denoted as the z axis. A uni- 
form THz radiation field (RF) E 1Ez (t) = Ecos(fif) = 
(E x , E y , 0) cos(fM) is applied in the x-y plane with the pe- 
riod To = 2tt/Q. The angle between the electric field and 
x axis is #e- By using the Coulomb gauge, the vector and 
scalar potentials can be written as A(i) = — Esin(f}i)/f2 
and (j)(t) — 0, respectively. We assume the well width 



with = E^tf 1 + 75 6 f ), = £,( 7 «« - 3 75 e 4 
7c H = E z ^!f l - It i s noted from these equations that the 
magnitude of the Rashba term can be tuned by the exter- 
nal electric field E z applied on the sample. 15 ' 16 In Eqs. 
©-©, ilf h , i\f h , ifi\ ifi\ itf 1 are the Rashba 
coefficients. 13 They depend both on the property of ma- 
terial and QW well width. 

Similar to Refs. d and [l(| by employing the Floquet 
theorem^ the solution of the Schrodinger equation with 
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time-dependent Hamiltonian H\(K,t) can be written as 

3> A (k ^ _ e -i{(B^+B^„ l )t-&ok-E[cos(Ot)-l]-7sin(20t)} 

oo 

xe ^ A (k) t <Pl s (k)e innt (6) 
n=— oo 

Here s = ± represents the two helix spin branches; 
= k 2 /2m x is the kinetic energy of HHs or LHs; 



E em — e 2 E 2 / (4m* x il 2 ) is the energy induced by the RF 
due to the DFKE; b = e/(m* x n 2 ); 7 = E em /(2Sl). 

<.(k) = = Cl^J) in Eq- © are the ex - 

pansion coefficients of the Floquet states with a = 1 (—1) 
representing spin- up f (-down J.) in the laboratory coor- 
dinates (along the z axis). <k(k) is the corresponding 
eigenvalue and can be determined by 



[nfi - <7 s A (k)]<^ s + {[£& (k) ± iaD^(k)\ + 2(e/2n) 2 [D^(k) ± «r£&(k)]}0-* 
+ {i(e/2fi)[I>* (k) ± *a7J A 2 (k)] + 3i(e/2n) 3 [I&(k) ± *a# 3 A 2 (k)] } (^ 1>s - 

- (e/2fi) 2 [I>* (k) ± ivDUm^s + 4>n-2,s) ~ i(e/2fi) 8 [D*(k) ± lf x J D A 2 (k)](^ s - = , (7) 

where 

r>HH _ HH .3 , ~HH k 2 , D HH _ HH ,3 , HH j.2 . 

= 7 a L % 3 + lt H k% + ^ H (k 2 z )k y , Dg = rt H kl + rf»k 2 y k x + ^ H (k 2 z )k x , 

D" H = 3ja H k 2 y E y + ^ H (k 2 x E y + 2k x k y E x ) , D» H = 3j™k 2 x E x + 7 ™ \k 2 y E x + 2k y k x E y ) , 

D^ H = ila H k 2 y E y + jt H (k 2 x E y + 2k x k y E x ) + ^ H (k 2 z )E y , 7J A 2 = 3^k x E 2 x + ^{2k y E x E y + k x E 2 y ) , 
D A = 3^k y E 2 y + ri(2k x E y E x + k y El) , D*f = 3^ H k 2 x E x + ^ H (k 2 y E x + 2k y k x E y ) + ^ H (k 2 z )E x 

^3 A i = 7a A ^ 3 + lbE 2 x E y , D$ 2 = >$E* + ^E 2 y E x . 



All eigenvalues can be written as q s . n = Is. o + nO. where 
g S; o is the eigenvalue in the region (— 0/2, Q/2]. It is 
evident that q SyTl and g Sj o are physically equivalent. We 
also find s = + branch and s = — branch satisfying the 
relations: 



g_(k) = -«+(k) 



(8) 
(9) 



With the help of Green function, we can calculate the 
density of states (DOS) p a ^ a and the density of spin po- 
larization(DOSP) p„- a ^- 
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S=± H<h — 



s i(n-m)fiT 



xi? criiCr2 (s;n,TO;k)J; 1 (-26 k • Esin(OT)) 
x J /2 (2 7 cos(2r!T))<5(w - [E k + E em 
-(li+2f a +n + m)fi/2 + g,(k)]) , (10) 

in which J n (x) is the Bessel function of nth order, 

Rai,a 2 (s;n,m;k) = fao-j Vn,a(k))(0m,s' (k)f?<r 2 ) with rj a 
standing for the eigenfunction of a z . It is seen from Eq. 
(TTU|) that these densities are periodic functions of T with 
period T . The DOSP is nonzero only when both the 
RF and the SOC are present. Furthermore, the induced 



magnetic moment can be written as 



M(T) 



(M x (T),M y (T),M z (T) 
2gp B r^m 



n T + n L J_ QO 
^(Phl ~ Pl,l)) 



dw(Rep Ta ,-Imp T 



1- 



(11) 



where the Fermi energy Ep(T) is determined by n a — 
J^ F p^c^jT) du where n a represents the hole con- 
centration. Eq. (TTTj) has been simplified by using the 
fact that p ai ,a 2 — P% 2 <y\- 1^ ^ s evident that Ep(T) and 
M(T) both oscillate with the period To- Due to time re- 
versal symmetry, the DOSP is an odd function of the 
time pa- t -a-(T,Lu) = — p a ,-<j{— T, lu), and therefore the 
DOSP averaged over time reduces to zero. Besides, the 
DOS is an even function p a _ a (T,u)) = p ata (—T,uj) and 
Pl^(T,u>) = p^i(T,uj), thus the RF in the x-y plane 
cannot induce magnetic moment along the z axis. These 
characters are similar to those of a 2DEGi2 

We numerically solve the eigen-equation Eq. ([7]) and 
calculate the DOS and the DOSP through Eq. flU]). 
One can further obtain the magnetic moment by us- 
ing Eq. (fTTj) . In the calculation we choose a = 10 nm, 
E z = 30 kV/cm. The material parameters of GaAs are 
as follows)! 3 .^ 7i = 6.85, 72 = 2.1, 73 = 2.9, A = 0.341 



eV, g LH = 1.2, g H H = 3.6, m 



LH 



0.0537 mo and 
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m* HH — 0.171 mo. In order to ensure the validity of 
the model which we adopt, we must keep the highest 
sideband of HH well separated from the lowest sideband 
of LH. Moreover, the HH and LH bands can be splittcd 
by 50 meV by adjusting the applied strain^ According 
to these, we choose E = 0.1 kV/cm, f2 = 0.1 THz in the 
following calculation. 
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FIG. 1: Time-averaged DOS of HHs under THz field with 
E = 0.1 kV/cm, Q = 0.1 THz (solid curve) and without THz 
field (dashed curve). 
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FIG. 2: DOSP of HHs for (a) 6 E = and (b) E = tt/4 at 
T = T / 4 with E = 0.1 kV/cm, fi = 0.1 THz. The real part 
and imaginary part are plotted as dashed curves and dotted 
curves respectively. The solid curves is Y as a function of Ef- 
Note the scale of Ef and Y are on the upper and right frame 
of the figure. 
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FIG. 3: Magnetic moment M of HHs versus time for 0b = 0, 
E = tt/6 and E = tt/4 with E = 0.1 kV/cm, tl = 0.1 THz. 
The concentration of HH is n-\ = nj. = 0.7 x 10 11 cm" 2 and 
Ef is about 0.35 meV. 



In Fig. [T] we compare the time-averaged DOS with and 
without the THz held. Due to DFKE, the main absorp- 
tion edge has a blueshift and the DOS becomes finite 
below the band gap. The DOSP at T = T Q /A are plot- 
ted in Fig. H] for THz field along two different directions, 
(a) 6 E = and (b) E — n/4. It is seen from Fig. [2(a) 
that only the imaginary part of DOSP is finite. From 
Eq. I|lip. we can find that the induced magnetic mo- 
ment is along the y axis. This is similar to the 2DEG 
case with Rashba SOCj£ In Fig. [2(b), we can see that 
Rc/jj^ = — Impf.j. Thus the induced magnetic moment 
is along the (l,l,0)/v2 direction, i.e., parallel to the 
THz field. These results indicate that the direction of 
the induced magnetic moment varies with that of the 
THz field, which is different from the 2DEG case with 
Rashba SOC. 9 This is due to the anisotropy of the SOC 
Hamiltionian. 

In Fig. [3l the magnetic moment of M is plotted as 
function of time for Oe = 0, Oe = 7r/6 and 0e = tt/4 
with E = 0.1 kV/cm, Q, = 0.1 THz. It is noted that 
that the magnitude of magnetic moment depends on the 
direction of the THz field 6e- The magnetic moment is 
the smallest for 0e = and the largest for 0e = vr/4. 

In Fig. [2J we also plotted T as function of Ep , where 
T = /** Im p Ti X cbj. It is noted that T is very small for 
large enough Ep, hence M becomes negligible when the 
concentration of the hole gas is high. This can be under- 
stood as follows: By interchanging the order of integral, 
one has 

/oo 
cik Jh (-26 k ■ E sin(fiT)) 

x ]T # u ( S ;n,m;k)e^- m )" T 

s,m,n 

x ^ J /2 (2 7 cos(2fir)) . 
h 

By virtue of Eqs. © and ©, one gets R^^ (s; n, m; k) = 
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—R^l(—s; —to, — n; k). Thus the terms of s = + branch 
compensate those of s = — branch, and the integral of 
the DOSP over the whole range (—00, 00) is zero. On the 
other hand, the DOSP decays to very small value with 
increasing u due to the interference caused by the mo- 
mentum dependence of the SOC. Hence the contribution 
to the magnetic moment at large u) is negligible. Accord- 
ingly, M becomes very small when Ep is large, i.e., the 
hole concentration is high. Our calculation shows that 
this is also true for 2DEG with Rashba SOC. 

In conclusion, we study the effects of the intense THz 
field on 2DHG with Rashba SOC. We calculate the DOS 
and DOSP. We also show that the a THz magnetic mo- 
ment can be excited for low hole concentration. It is 
noted that the direction of the THz field has a strong 
influence on the angle between the induced magnetic mo- 
ment and the THz field, as well as on the amplitude 



of the magnetic moment, which is quite different from 
2DEG with Rashba SOC case. We also point out that 
the magnetic moment becomes very small if the hole con- 
centration is high enough, due to the interference caused 
by the momentum dependence of the SOC. This effect 
also appears in 2DEG. 
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